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A relativistic transport model is used to study H“ production in 2 — 11A GeV Au+Au col¬ 
lisions. Introducing the strangeness-exchange reactions between antikaons and hyperons as the 
sources for E“, we find that the cascade yield in these collisions is in reasonable agreement with 
the data. Although the E _ abundance does not reach chemical equilibrium unless the cross section 
for strangeness-exchange reactions is enhanced by six times, it exhibits the strongest enhancement 
with increasing centrality of collision and with increasing beam energy. 
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I. INTRODUCTION 


One of the most challenging problems in relativis¬ 
tic heavy-ion collisions is to explain the observed en¬ 
hancement of singly strange and multistrange hadrons. 
The shorter equilibration times for strange quarks in a 
quark-gluon plasma (QGP) than the strange hadrons in 
hadronic matter had led to the conjecture [1] that en¬ 
hanced production of (multi)strange hadrons might be 
an indication of QGP formation in the early stage of 
the collision. However, a firm conclusion of strangeness 
enhancement as a signal for the QGP can only be es¬ 
tablished if baseline hadronic processes fail to explain 
the strange particle abundances. For the singly-strange 
kaons and lambdas, it was found that the observed en¬ 
hancement at both the AGS [2] and the SPS [3] ener¬ 
gies could indeed be understood within a hadronic sce¬ 
nario [4]. This scenario, on the other hand, cannot 
explain the enhanced production of the multistrange 2 
and Yl (and their antiparticles) observed in Pb+Pb col¬ 
lisions at SPS energies [5,6]. The latter enhancement 
has thus been associated with production mechanisms 
which include possible QGP formation [7]. The doubly 
strange 2 has also been measured in heavy ion collisions 
at the AGS [8]. Compared to that expected from pri¬ 
mary nucleon-nucleon collisions, the observed enhance¬ 
ment is even more pronounced than in the SPS experi¬ 
ments. Although this is not surprising as experimental 
data have shown that the ratio of strange to nonstrange 
particles (such as K + /tt + ) in A+A relative to p+p col¬ 
lisions increases gradually from RHIC energy to lower 
beam energies at SIS [9], it is interesting to find out if the 
measured 2 yield in this collision is consistent with the 
hadronic or the QGP scenario. In this paper, we shall 
study the production of multistrange 2“ from Au+Au 
collisions at the AGS energy of 6 A GeV using the rel¬ 
ativistic hadronic transport model ART [10]. We shall 
show that strangeness-exchange reactions between the 
antikaon and lambda(sigma) in hadronic matter alone 
can explain the measured 5 yield [8], and exotic pro¬ 
cesses for multistrange particle production at the AGS 


energies are thus not required. 

II. THE ART MODEL AND THE 
STRANGENESS-EXCHANGE REACTIONS 

The ART model is a hadronic transport model that in¬ 
cludes baryons such as N, A(1232), W(1440), iV*(1535), 
A, E, and mesons such as 7r, p, lo, rj, K, K*. Both 
elastic and inelastic collisions among most of these parti¬ 
cles are included by using the experimental data from 
hadron-hadron collisions. The ART model has been 
quite successful in explaining many experimental obser¬ 
vations, including the surprisingly large kaon antiflow 
[11,12] in heavy ion collisions at AGS energies. The 
ART model also allows us to understand whether or 
not the strongly interacting matter formed in these col¬ 
lisions reaches chemical and/or thermal equilibrium. In 
the present study, we extend the ART model to include 
perturbatively the 2 particle as in the studies for other 
rare particles using the transport model [6,13,14]. 

In hadronic matter, the 2 particle is mainly produced 
via the strangeness-exchange reactions Jv(A, S) —► 7 t 2 
because the contribution from the associated produc¬ 
tion reactions NN —> ENKK and 7 t(A, E) —> KE 
are Okubo-Zweig-Iizuka suppressed. Since there is no 
empirical information on the strangeness-exchange reac¬ 
tion for 2 production, we use the cross section obtained 
from the gauged flavor SU (3)-invariant Lagrangian in the 
coupled-channel approach [15]. Because the large thresh¬ 
old for final states with the rj particle, only the reactions 
K( A, E) —+ 7t 2 are considered here. Writing the spin- 
and isospin-averaged cross sections for these reactions as 

a KA^irE = T~|-^A'A^-n-s| 2 ! 

4 PA 

a KY.^irE = ~Tcf 7 l-^AS—>7rH P > ( 4 ) 

1Z p K 

where pj^ and p are initial antikaon and final pion mo¬ 
menta in the center-of-mass system, the theoretical cross 
sections of Ref. [15] are found to be well fitted by the 
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following squared invariant matrix elements 
and \M R x^ s \ 2 [16]: 

\M Ra ^ s \ 2 = 34.7^ mb, 

|-^iC£—>7rH| I" = 318 (l “) 0 6 (tT") 17 mb - ( 2 ) 

1 /2 

In the above, the threshold energy s 0 in the center-of 
mass system is 1.611 GeV and 1.688 GeV for the reac¬ 
tions KA —> 7 rS and KT, —► 7rS, respectively. We note 
that except near threshold, these cross sections are of the 
order of 5-10 mb. 

The cross sections for the inverse reactions 7rS —> KA 
and 7 rS —> K E, which are needed for treating 3 anni¬ 
hilation, are related to those for 3 production by the 
principle of detailed balance, i.e., 

- 1*4 

*+rH— >K A o 2 

A 7+ 

^TrE-tiCE — —1 + 

pi 

We have also included the production of the resonance 
5* (1530) via the reaction 7 t 5 —» 3* (1530), and its cross 
section is given by a Breit-Wigner form with the empir¬ 
ical width of T s . = 9.5 MeV. The cross section for the 
inverse reaction of 3*(1530) decay, i.e., 3*(1530) —> 7 t3 , 
is also obtained from the detailed balance relation. 


III. RESULTS 


A. Time evolution of hadron abundances 



FIG. 1. Time evolution of central density and hadrons at 
midrapidity |y| < 0.5 in 6A GeV Au+Au collisions at impact 
parameter b < 3 fm in the ART model. 


Figure 1 displays the time evolution of midrapidity 
hadrons and the central density obtained in the ART 
model for central Au+Au collisions at 6 A GeV. The 
abundances for baryon resonances (left panel) attain 
peak values at the maximum compression stage of col¬ 
lision at a time t « 6 fm/c after the initial contact 
of the colliding nuclei. During subsequent decompres¬ 
sion, the resonances decay to produce the stable nucle¬ 
ons and pions. The kaons produced from baryon-baryon 
and pion-baryon collisions in ART model, tends to sat¬ 
urate soon after the maximum compression stage (right 
panel). On the other hand, the antikaon yield steadily 
decreases during the expansion stage due to their strong 
absorption with the abundant nucleons at the AGS en¬ 
ergies via KN —> 7t(A, E). The hyperons (A, E,S plus 
their respective resonances) are produced within 8 fm/c 
when the density is high. Thereafter, the abundances 
are rather insensitive to hadronic scatterings, if any. We 
find that the strangeness-exchange reactions between an- 
tikaons and hyperons lead to an appreciable production 
of 3. 




FIG. 2. Time and transverse radius distributions of midra¬ 
pidity (|t/| < 0.5) hadrons at freeze-out in 6A GeV Au+Au 
collisions at impact parameter b < 3 fm from ART model. 

The time and spatial distribution of midrapidity par¬ 
ticles at freeze-out are shown in Fig. 2. Most of the S’s 
freeze-out between 3 and 7 fm/c, whereas the freeze-out 
distribution of protons and pions is centered at about 12 
fm/c (left panel). Due to the fact that strange hadrons 
have much smaller scattering cross sections with nucle¬ 
ons and pions, they decouple quite early from the sys¬ 
tem. The dominantly soft collisions among the non- 
strange hadrons lead to a long tail in the nucleon and 
pion distributions. The transverse radius distribution 
(right panel) indicates that the 3 and other hyperons 
have source sizes similar to the initial source. 

B. Impact parameter dependence 

In Fig. 3, we show the impact parameter dependence 
of the total yield of different particle species obtained 
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from the ART model for Au+Au collisions at 6 A GeV. 
Except for very peripheral collisions, a K~/K + ratio of 
0.09 is obtained and found to be independent of colli¬ 
sion centrality. This trend is also observed experimen¬ 
tally from ~ 11A GeV up to RHIC energies [17,18]. 
Within a statistical approach, this implies that K + and 
K~ have similar freeze-out volumes. Indeed, in the 
transport calculation kaons and antikaons exhibit sim¬ 
ilar mean freeze-out times (~ 8.3 fm/c) and transverse 
radii (~ 4.5 fm) for central collisions. The A + E° ex¬ 
hibits a similar impact parameter dependence because 
of its associated production with kaons via the reaction 
( 7 r, p, u>){N 1 A, N*) —> K( A, E). Except for the very cen¬ 
tral collisions the A + E° yield in the ART model is found 
to be consistent with the E895 data. 



b (fm) 


FIG. 3. Calculated impact parameter dependence for the 
total yield of K + , K ~, A + E° (top panel) and (bottom 
panel) in Au+Au collisions at 6A GeV. The dependence of 
E~ production on the cross section a for I\ (A, E) —* nE is 
also shown. The E895 data for the A + E° (top panel) and S 
yield (bottom panel) are shown by solid circles. 


The 2“ yield obtained by strangeness-exchange reac¬ 
tions between antikaons and lambda/sigma for different 
impact parameters is also shown in Fig. 3. Similar to 
the singly strange hadrons, the 5“ yield exhibits nearly 
a linear dependence on impact parameter. This is due to 
the fact that 2“s are mostly produced from secondary 
collisions in the ART model. Fig. 3 further shows that 
while the A+E° is enhanced by a factor of 4, the 2“ yield 
grows steadily by a factor of 6 from peripheral to central 
collisions. This specific hierarchy of strangeness enhance¬ 
ment of 5 > A+E° has been reported for Pb+Pb collision 
at the SPS by the WA97 collaboration [19]. However, at 
SPS energies, the 2“ yield tends to saturate for a large 
number of participants V part > 100, and its enhancement 
(by ~ 3) in central collisions is smaller than our predic¬ 
tion for AGS energies. The steady growth of 2“ with 


centrality is a clear indication that its yield is sensitive 
to the volume of the fireball formed in the collision. Fur¬ 
thermore, we find that the 2“ yield from strangeness- 
exchange reactions in the ART model is in reasonable 
agreement with the E895 data [8]. This is in contrast 
to SPS energies, where the cascade from the transport 
model was found [6] to underpredict somewhat the WA97 
data, especially for central collisions. 


C. Comparisons with the statistical model 


It may be interesting to find if strange hadrons also 
reach chemical equilibrium in heavy ion collisions at 6 A 
GeV. To this end we compare results from the transport 
model with those from the statistical model based on the 
grand canonical ensemble [20,21] with complete thermal, 
chemical, and strangeness equilibrium. In this model, the 
particle density in the Boltzmann approximation is given 
by 


Ni 

V 


grnjT 
27r 2 


K 2 (mi/T)e^ T , 


( 4 ) 


in the usual notation [20,21]. Assuming 2 to be also in 
chemical equilibrium, the cascade to nucleon ratio is 


Ne = (ms \ 2 K 2 (mz/T) , 

N n \m N J K 2 (m N /T) 

- i ( Na V ( msmN \ 2 K 2(m~/T)K 2 (m N /T) 

\N n ) V m \ J [K 2 (m\/T )] 2 

For central Au+Au collisions, we use the ART model 
result for the ratio N\/Nn = 0.025. This gives a 2~ 
yield of N s - = 0.35 — 0.38 at a freeze-out temperature 
of T = 120 — 170 MeV. Note that the production of 
the strange particle 2 is less than 1 per event. Thus, 
strangeness conservation must be implemented locally, 
i.e., a canonical instead of grand canonical ensemble must 
be used [22]. Therefore the present estimate of 2“ gives 
only an upper bound. The higher value for the estimate 
of 2“ from the statistical model compared to that ob¬ 
tained in the present transport approach indicates that 
2~ does not reach chemical equilibrium. Furthermore, 
like other strange hadrons, 2 freezes out relatively earlier 
than nonstrange hadrons with a mean time and trans¬ 
verse radius of 8.4 fm/c and 4.6 fm when the central 
energy density is about e ~ 1 GeV/fm 3 . 

The 2“/(A + E°) ratio of 0.026 for b = 0 — 3 fm in the 
ART model is smaller than the equilibrium model esti¬ 
mate of 0.042 but slightly larger than the E895 data [8]. 
At midrapidity, the 2“/(A + E°) ratio of 0.12 has been 
observed by the E810 Collaboration for central Si+Au 
collisions at 14.6A GeV [23]. Since 2 is generated from 
A, E by strangeness-exchange reactions, its yield would 
naturally be more peaked at midrapidity than that for 
As. This is corroborated by a higher 2“/(A + E°) ratio 
of 0.03 at midrapidity \y\ < 0.50. 
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Another way to see if the yield of 5 approaches that 
for chemical equilibrium is to study how it changes when 
the cross section is artificially increased for the reaction 
K( A, E) ^ tt 5. This method was introduced in Ref. [24] 
to study whether or not K~ reaches chemical equilibrium 
at energies available at SIS/GSI. The latter energies are 
below the threshold for K~ production from nucleon- 
nucleon collisions in free space. Figure 3 shows that the 
5“ yield increases rather dramatically with increasing 
cross section a. Equally important is the fact that for a 
given o', the enhancement is even more pronounced for 
central collisions. Nonetheless, it is only for very large 
cross sections (~ 6er), where the hydrodynamical limit 
may have been reached, that the 5“ approaches chemi¬ 
cal equilibrium with A+- = 0.55. 

Interestingly, this value of 5~ is even larger than that 
predicted from the statistical model. Possible reasons for 
this result are: (i) The temperature at complete chemical 
equilibrium, if reached at all, is much higher than that 
in the statistical model; (ii) The increased cross section 
has altered the reaction dynamics especially for K and 
A, i.e., a large abundance of Ks and A + E°s so produced 
will increase the 5 at a higher temperature; or (iii) large 
numbers of 5s, produced at an early and dense stage do 
not have sufficient time to be annihilated by rescattering 
because of rapid expansion of the system during its sub¬ 
sequent evolution. The possibility of large 5 production 
from case (i) alone is remote as its yield in the present 
statistical model is found to be rather insensitive to tem¬ 
perature. 


D. Excitation functions 



FIG. 4. Excitation function for the total yields of K + , 
A' - , A + E° and H~ (top panel) and the ratios I\~ / K + , 
E _ /(A + E°) and S~/K~ (bottom panel) in central (6 = 0 — 3 
fm) Au+Au collisions from ART model. 


In Fig. 4, we show the excitation function for 
the strange particle abundances (top panel) in central 
Au+Au collisions from the ART model. As expected, 
with increasing beam energy the multiplicities of the 
hadrons increases. The enhancement is more pronounced 
for antikaons and cascades as the energy available for 
these massive and rare particle production increases. 
This is clearly evident in the bottom panel of the Fig. 
4 where the I\~ / K + ratio is found to increase monoton- 
ically with beam energy. The 5“ has the largest rate 
of increase with energy resulting in an increasing of the 
ratios 5 _ /(A + E°) and 'Br/K~ with beam energy. 


IV. SUMMARY 

In summary, we have used a relativistic hadronic trans¬ 
port model to investigate the doubly strange cascade pro¬ 
duction. We find that the strangeness-exchange reactions 
between antikaons and hyperons lead to substantial pro¬ 
duction of 5“ in Au+Au collisions at 6 A GeV that is 
in reasonable agreement with the data. This is in stark 
contrast to that at SPS energies, where the strangeness- 
exchange reactions underpredict the data for the cas¬ 
cade yield in central collisions. The success of the purely 
hadronic scenario eliminates the scope of an exotica that 
may be necessary to explain the 5 yield at the AGS en¬ 
ergies. Among all the strange hadrons, the 5 particle 
abundances reveal the strongest enhancement with cen¬ 
trality of collision and also with increasing beam energy. 
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